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Abstract 

The role of restorative coupling on synchronization of coupled identical harmonic oscillators is 
studied. Necessary and sufficient conditions, under which the individual systems’ solutions converge to 
a common trajectory, are presented. Through simple physical examples, the meaning and limitations 
of the theorems are expounded. Also, to demonstrate their versatility, the results are extended to 
cover LTI passive electrical networks. One of the extensions generalizes the well-known link between 
the asymptotic stability of the synchronization subspace and the second smallest eigenvalue of the 
Laplacian matrix. 


1 Introduction 


Studying the collective behavior of coupled harmonic oscillators has been a rewarding enterprize for 
researchers who try to enhance their understanding on a much-encountered phenomenon in nature: syn¬ 
chronization. For instance, it has been observed that two or more identical pendulum^] connected by 
means of dampers eventually swing in unison even if initially they are not synchronized; see Fig. [l] This 
outcome is not difficult to reach by intuition. Since the energy of the system can only leak out through 
the dampers, the pendulums should eventually settle to a constant energy state where there is no leakage. 
No leakage implies that the relative velocities are all zero. In other words, all the pendulums are moving 
at equal velocities at all times. This is only possible when they are synchronized. 



Figure 1: Damper-coupled pendulums. 


The simple example above has served as a starting point for many significant generalizations. In [5] 
Ren studies synchronization of coupled harmonic oscillators allowing time-varying oscillator dynamics as 
well as time-varying and asymmetrical dampers. The case where the damping between a pair of oscillators 
becomes effective only when the two are close enough is investigated in [10]. The effect of nonlinear 
damping is analyzed in j3] and of impulsive damping in [14]. A sampled-data approach is adopted in 
mm- Adaptive damping is covered in [S] and synchronization in the presence of noisy damping is 
considered in m- Note that all these works consider only dissipative coupling (e.g. dampers). From 
the engineering point of view this choice is not surprising because introducing restorative coupling (e.g. 
springs) will in general deteriorate performance by causing longer and more oscillatory transient behavior; 
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1 We restrict our attention to the small oscillations, where the pendulum can be represented by a linear model. 


1 














for instance, the simulation results show that the three pendulums in Fig. [2] synchronize much less rapidly 
than those in Fig. [1] Perhaps this may partly explain why collective behavior of spring-coupled oscillators 
has attracted more physicists than engineers. While for the engineer spring is an option to couple two 
units, for the physicist it represents an inherent characteristic of interaction. Relevant investigations 
in the physics community goes as far back, if not further, as the work of Fermi et al. [4] where chain 
of nonlinearly coupled oscillator-like particles were studied. Due to the richness of the subject and the 
increasing variety of applications in both inanimate and biological systems, the area has maintained its 
livelihood throughout many decades. See, for instance, 0, h, © 0 for recent progress. 



Figure 2: Damper- and spring-coupled pendulums. 


Through this paper we aim to provide a comprehensive analysis of the collective behavior of identical 
harmonic oscillators coupled by both restorative and dissipative components. To the best of our knowledge 
a detailed treatment of synchronization has not yet been reported for this setting, where two different 
interconnection graphs are simultaneously at work: the graph representing restorative coupling and the 
graph representing dissipative coupling. We present a necessary and sufficient condition on the associated 
pair of Laplacian matrices, under which the individual systems tend to oscillate in unison. We also point 
out a certain sufficient-only, yet easier-to-check set of conditions guaranteeing synchronization and exercise 
them on some simple real-world examples for clarity. Later, we attempt to extend our approach to the 
analysis of linear electrical networks of identical oscillators (of arbitrary order) coupled through passive 
impedances. For such networks we establish a link between synchronization and the eigenvalues of the 
(complex) node admittance matrix. This seems to be a natural extension of the well-known connectivity 
condition in terms of the second smallest eigenvalue of the (real-valued) Laplacian matrix. 

2 Coupled harmonic oscillators 

Consider the array of q coupled harmonic oscillators 

9 9 

Zi -f- uj 0 Zi T ^ ( dij^Zi zj ) T ^ ( r'-j.'y (Zj Zj ) 0, i 1,2,...,^ (1) 

9 = 1 9=1 

where Zi £ R and ujq > 0 is the frequency of uncoupled oscillations. The symmetric weights dij = dji > 0 
and rij = rji > 0 respectively represent the dissipative and restorative coupling between the ith and 
jth oscillators. Note that without symmetry, i.e., either dij ^ dji or ^ rji, the solutions are not 
guaranteed to be bounded unless some extra assumption is made. We take da = 0 and ru = 0. In 
this section and next we search for conditions on the triple (wo, {dij}, (r,; 7 }) under which the harmonic 
oscillators © synchronize, i.e., | Zi(t) — Zj(t) | —>• 0 as t — >• oo for all i, j and all initial conditions. 

Let D, R € R 9 * 9 denote the weighted Laplacian matrices associated to the topologies described by 
the dissipative coupling {dij} and the restorative coupling {r^}, respectively. That is, 


' dij 

— di2 

d\q 


' £9 

-ri2 

■ ~ri q 

— ^21 

dij • 

~d2q 

, R = 

-r 2 1 

Ej r u ■ 

■ -r 2q 

^<?1 

— dq2 

dqj . 


. -^91 

-r q2 ■ 

Yjj r qj . 


Note that these matrices are symmetric positive semidefinite since dij = dji > 0 and = rji > 0. 

In particular, we can write z T Dz = ^2j >i dij{zi — Zj ) 2 and z T Rz = ’Yhj > i r ij{ z i z j) 2 > where 2 = 
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[z\ Z 2 • • • z,f £ I’. Let us now rewrite (0 as 


z + UqZ + Dz + Rz = 0. 

This, using x = [ z T z T ] T £ R 2q , allows us to obtain 


x = 


0 I q 

-{cu$I q + R) -D 


x =: <&x 


( 2 ) 


where I q £ R qxq is the identity matrix. Employing the symmetric positive definite matrix 



u}q I q + R 0 

0 I q 


we can establish the following Lyapunov equality 


$ T P + P$ 


0 0 
0 D 


Since the righthand side is negative semidefinite, each solution x(t) of the system J2]) is bounded. More¬ 
over, by Krasovskii-LaSalle principle, x{t) should converge to the largest invariant region contained in 
the intersection T> fl {a: : x T Px < x(0) T Pi(0)} where 


V ■- 



0 0 
0 D 



It turns out that the condition 


null 


R — XI q 
D 


C range l q for all X £ C 


(3) 


(where l q £ M 9 is the vector of all ones) guarantees that this largest invariant region is contained in the 
synchronization subspace 


S := range 


l q 0 

0 lg 


In other words: 


Lemma 1 Let ([3]) hold. Then and only then 


x{t) £ V for all t ==>• x(t) £ S for all t 


(4) 


where x{t) is the solution of the system (j2|) . 


Proof. We first establish © => Q. Let x(t) = [z(t) T z{t) T ] T be a solution of the system @ that 
identically belongs to V. This means Dz(t) = 0. Also, 



0 I q ' 


' 0 

0 


0 I q ' 

X = 

-{u%I q + R) 0 

x — 

0 

D 

X = 

— (Wg/g + R) 0 


which implies 


z + + Rf)z — 0 . 


(5) 


Let Ai, A 2 , ..., \ p be the distinct ( p < q) eigenvalues of R. Since R is symmetric positive semidefinite, 
these eigenvalues are real and nonnegative. Consequently, the matrix [ui^Iq + R] is symmetric positive 
definite with eigenvalues Wq + Ai, Wq + A 2 , ..., Wg + A p . Therefore (0) implies that the solution has the 
form 0 §23] 


p 

z(t) = ReJ2e jUkt fk 

k= 1 


( 6 ) 
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where u>k = -\/u)q + A^ are distinct and positive, and each £*, £ C 9 (some of which may be zero) satisfies 

0 = ([wg Iq + R] — LU k Iq)l;k 

= (R - a k I q )tk . (7) 

Note that the condition Dz(t ) = 0 and © imply 


D&= 0 


since are distinct and nonzero. Combining 0 and © we can write 


£fc £ null 


R — Xklq 
D 


( 8 ) 

(9) 


Suppose now © holds. Then © implies £/. £ range l g for all k. By © this readily yields z(t) £ range l q 
for all t. Consequently, z(t) £ range l q and x(t) £ S for all t. 

Now we show the other direction © ==> ©. Suppose that condition © is not true. Then we can find 
an eigenvalue A» > 0 of R and a nonzero vector £* £ R 9 satisfying £* ^ range l q such that D£* = 0 and 
(R — A */ g )£* = 0. Let w* = y/uiQ + A*. Using the pair (w*, £*) let us construct the function z* : R — »■ R 9 
as z*(t) = Re(e J “’*£*). This function satisfies the following properties. First, since £* ^ range 1 9 , we 
have 


z*(0) = £* $. range l q . (10) 

Second, since _D£* = 0, we have at all times 

Di,(i) = Re(jcv*e? u * t = 0 . (11) 

Third, since [R + (uq — w*)/ 9 ]£* = 0, we can write at all times 

z*{t) + (wg/ g + R)z*(t) = -w*2*(t) + (wg Iq + R)z*{t) = Re(e 3 ’ ta ’* t [.R + (wg -w*)/g]£*) = 0 
which together with m leads to 

z* (t) + Dz * (f) + (wg/q + Rk)z* (f) = 0 . (12) 

Let a;*(i) = [,z*(f) T i*(t) T ] T . It follows from (fl2l) that x*(t) satisfies © and hence is a solution of the 
system. By CD we can assert that the solution x*(t) belongs identically to V, but (flOl) tells us that x*{t) 
does not identically belong to S. That is, the condition © fails. ■ 


To the question asked at the beginning we can now give the answer: 

Theorem 1 The harmonic oscillators © synchronize if and only if © holds. 

Remark 1 Note that the condition © does not depend on the natural frequency uio ■ 

Before we end this section we attempt to interpret condition ©. Recall that, given matrices C £ R mx " 
and A £ R™ xn , the unobservable subspace of the pair (C, A) is 


unobs (C, A) = null 


G 

CA 

CA "- 1 


The below result reveals the meaning of condition © from the observability point of view. 

Theorem 2 Condition © holds if and only if unobs (D, R) = range l q . 

Proof. Suppose © fails. Then we can find an eigenvector £ of R satisfying £ ^ range l q and D£ = 0. 
Let A be the corresponding eigenvalue, i.e., i?£ = A£. We can write 


D 


to 


Q 

DR 

£ = 

DRf, 

= 

A Z?£ 

_ DR q ~ 1 


_ Di? 9_1 £ _ 


_ A 9-1 L>£ _ 
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Therefore £ £ unobs ( D , i?). Since £ ^ range l g , we must have unobs ( D , i?) ^ range l g . 

Now we show the other direction. Suppose unobs (D, i?) / range l g . By definition we have Dl q = 0 
and Rl q = 0, meaning unobs (D, R) D range l q . Consequently, 2 < dim unobs (D, R) =: t. Let 
{? 7 i, 772 , •••, r]{} with 771 = 1 ? be an orthogonal basis for unobs (D, R). Since R is symmetric its 
eigenvectors form an orthogonal basis for C 9 . Let this basis be {£ 1 , £ 2 , £ g } with £1 = l g . Note 

that we have span{ 772 , 77 ^} C span{£ 2 , £ g }. Now, let us choose an arbitrary nonzero vector 

w £ span { 772 , ..., 77 ^}. Since unobs (D, R) is R-invariant we have Rw £ span { 771 , ..., 77 ^}. Moreover, 
w £ span {£ 2 , ..., £ g } implies Rw £ span {£ 2 , ..., £ g } because £i are eigenvectors. Hence we can write 

Rw £ span { 771 , ..., 77 ^} n span {£ 2 , ..., £ g } 

= span { 772 , ..., ru} . 

This implies (since w was arbitrary) that span { 772 , ..., 77 ^} is R-invariant. Consequently, span { 772 , ..., 77 ^} 
contains at least one eigenvector £ of R. It must be that £ ^ range l g because l q ^ span { 772 , ..., 77 ^}. 
Let A be the corresponding eigenvalue, i.e., Rf = A£. Since £ £ unobs (D, R) we have £ £ nullD yielding 
null D n null (I? — A I q ) D span{£} which implies that J3]) fails to hold. ■ 


3 Sufficient conditions for synchronization 


Although condition ([3]) tells us definitely whether a given array of harmonic oscillators will synchronize 
or not, it may nevertheless be expensive or simply impossible to employ when, for instance, the number 
of oscillators is large or certain parameter values are unknown. Therefore it is worthwhile to look for 
sufficient-only, yet simpler-to-check conditions to determine synchronization. This is what we intend to 
do in this section. 

Recall that an undirected graph is a pair of sets (V, £) where V = (ci, 772 , •••, v q } is the set of 
vertices and the elements of the (possibly empty) set £ are some (unordered) pairs of vertices ( Vi , Vj). 
Let us now introduce two graphs associated to the array of harmonic oscillators © as follows. The graph 
Td = (V, £d) describes the interconnection associated to dissipative coupling and is such that ( Vi , Vj) £ £d 
when dij ^ 0. Similarly, T r = (V, £f) denotes the restorative coupling topology and (t>j, Vj ) £ £ T when 
Tij ^ 0. 

By construction Dl q = 0, which yields nullIA D range l q . Since a graph is connected when the 
eigenvalue of the associated Laplacian matrix at the origin is simple, we have nullD = range l q when Td 
is connected. Note that © is trivially satisfied if null D = range l q . Therefore we can assert: 

Corollary 1 The harmonic oscillators © synchronize if the dissipative coupling graph Td is connected. 


As mentioned earlier, a collection of identical pendulums connected (only) by dampers eventually 
synchronize. What the above result adds to this statement is that even if we supplement the collection 
by springs connecting some pairs of pendulums, the tendency for synchronization cannot be destroyed. 
Corollary [T] is hardly surprising. Now we move on to establishing a less evident result. We begin by 
defining the matrix Ra £ M 9X? as 



' Ej tv 

_fi2 

■ -h q 


-T21 

Ej j ■ 

■ ~r 2q 

Ra = 






~fq 2 

Ej r <n 

We let Ts = 

(V, &) 

where £s = 

£r U £d- 


where 


_ J ' y 


for d^ = 0 , 


0 for d^ ^ 0 . 


Likewise, employing the set difference £a = £ r — £d we 
define the graph Ta = (V, £a)- Let Ta have c connected components, which we denote by T^ = (Vf, £f) 
for i = 1, 2, ..., c. By definition all the pairs (Vf, 14) and (£e, £k) are disjoint for £ ^ k. Moreover, 
(J/14 = V and (J f £{ = £a- Note that if Ta itself is connected then c = 1 and Ti = Ta- Let nf, = |14| be 
the number of vertices that belong to T^. Without loss of generality let the vertices V{ be such labelled that 
Vi = {vi, U 2 , ..., 77 ni }, V 2 = {?7 ni+ i, 77 ni+2 , ..., 77„ 1+ri2 }, and so on. Then Ra has the block diagonal 
form 


Ra 


Ri 0 

0 f? 2 


0 

0 


0 0 R c 
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with Re £ M. neXne . Note that each Re is symmetric positive semidefinite and satisfies Reln e = 0 meaning 
there is an eigenvalue at the origin. Since is connected this eigenvalue at the origin is simple. Hence 
the eigenvalues of Re can be ordered as 0 = Xij < \ 2 ,e < ■ • ■ < A . Now, for each £ = 1, 2, ..., c, 
define the following system 

r)e + (wq I ne + Re)r]e = 0 (13) 

with ije = [zat+i z ai +i ''' Zae+n e ) T £ where o\ = 0 and ae = ni+r^d- \~ne~i for £>2. Note that 

we can write z = [z± Z 2 • • • z q ] T = [ryf rfi ■ ■ ■ vI) T ■ Let fie = {w > 0 : w 2 = uiQ + Xk,e , k = 1, 2, ..., ne } 
denote the set of characteristic frequencies of the system m- Note that the frequency wo belongs to 
every Vtg because X±,e = 0. Now we list a sufficient set of conditions guaranteeing synchronization. 

Assumption 1 (P) The harmonic oscillators (HD satisfy the following conditions. 

1. For all l — 1, 2 , ..., c the system (ITdl) is observable from each Zk for all k £ {i : Vi £ Ve}. 

2 . n 1 nn 2 n---nn c = {w 0 }. 

3. is connected. 

A more mathematical (less physical) version of Assumption |T| reads: 

Assumption 1 (M) The harmonic oscillators (HD satisfy the following conditions. 

1. None of the matrices R±, R 2 , ..., R c has an eigenvector with a zero entry. 

2. A = 0 is the only common eigenvalue of the matrices R±, R 2 , ..., R c . 


3. null R f~l null D = range l q . 

Theorem 3 The harmonic oscillators HD synchronize if Assumption [7] holds. 


Proof. If we can show that Assumption |T) implies condition ([3]) then by Theorem [T] the oscillators 
must synchronize. Let us establish the implication by contradiction. Suppose that © is not true but 
Assumption |T] holds. Then we can find an eigenvector £ ^ range l q satisfying D£ = 0 and (R — A I q )f = 0 
for some A. This eigenvalue A cannot be zero for then we have £ £ null I? D null D which contradicts the 
third condition of Assumption [T| Let us therefore study the case A 7 ^ 0 in the sequel. 

Let us employ the partitions £ = [z\ Z2 • • • z q ] T = [rff 77^ • • • ??J] T . Since Df = 0 we can write 
0 = £ T Df = '}2j >i dij(zi — Zj) 2 which implies that for a given pair (i, j) of indices either dij = 0 or 
Zi = Zj. Since dij = 0 means fij = rij we have fij{zi — Zj) = rij(zi — Zj) for all i, j. Now we can proceed 
to claim Rf = R/\f because 


" T,j r lj{ z l-Zj) ' 


' EjhM-Zj) ’ 

. E j r qj{z q -Zj) _ 


. J2j f qj(z q - Zj) _ 


= #a£- 


Therefore I?a£ = A£. Then the block diagonal structure of Ra allows us to write Rerje = Xpe for all 
t = 1, 2, ..., c. Since by assumption nonzero eigenvalue A is not common to all Re we should have rje = 0 
for at least one index t. Also, again by assumption, no entry of r\e can be zero whenever r\e 7 ^ 0. That 
is, if r\e 7 ^ 0 then Zk 7 ^ 0 for all k £ {i : Vi £ Ve}. Let us now define two (nonempty) sets of indices 
I := {i : Vi £ Ve, rje 7 ^ 0, £ = 1, 2, ..., c} and its complement J := {1, 2, ..., q} — T. Note that Zi 7 ^ 0 
for i £ I and Zi = 0 for i £ J. For any pair of indices (i, j ) with i £ I and j £ J we can assert the 
following, (i) dij = 0. Because Zi — Zj 7 ^ 0 and 0 = £ T Z?£ > dtj{zi — Zj) 2 . (ii) fij = 0. Because, by how 
we constructed the sets I and J , the vertices Uj and Vj cannot belong to the same vertex set Ve- Then by 
the block diagonal form of Ra the entry fjj must be zero, (iii) r,;. ; = 0. Because dij = 0 means = fij. 

Construct the vector £ = [z\ £2 ■ ■ ■ z q ] T with entries Zi = 0 for i £ T and Zi = 1 for i £ J . Clearly, 
£ ^ range l q . We can write 


£ r D£ = ^dijizi-zj ) 2 
j>i 

= — dij(zi-Zj) + dij(zi-Zj) +- ^2 dij(2i — Zj) 

= 2 5 Z “ °) 2 + 2 5 Z ~ !) 2 

i,jej 

= 0 
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where we used the fact that d, ;/ = 0 when jgl and j € 3 ■ Then £ T D£ = 0 implies D£ = 0 because D is 
symmetric positive semidefinite. Since we also have that r,j = 0 when i 6l and j G 3 , we can similarly 
establish = 0. Hence £ G nulli? (~l null I?. But £ ^ range l q . This contradicts the third condition of 
Assumption |T] ■ 

As stated earlier, Assumption[l]is only sufficient for synchronization. Hence if an assembly of harmonic 
oscillators fail to synchronize, at least one of the three conditions listed therein must not hold. Of those 
three conditions, the necessity of the third one (that Ts is connected) is evident. However, the relation 
of the remaining two conditions to synchronization is subtle and requires due attention. To better 
understand the meanings of those conditions we now provide two examples, where harmonic oscillators 
do not synchronize. Each example violates one of the first two conditions of Assumption [T] 


Example 1 Consider the following four coupled harmonic oscillators 

Z\ + UJo z l + r ( z i ~ z 2) = 0 

z 2 + Wq -2 + d(z 2 — z±) + r(z 2 — Zi) + r{z 2 — 23 ) = 0 
Z 3 + Wq-3 + r(z 3 - Z2) = 0 

Z4 + LJqZ4 + d(z4 — z 2 ) = 0 


where Wq, d, r > 0. The associated D and R matrices are provided below. 


' 0 

0 

0 

0 ' 


r 

—r 

0 

0 ' 

0 

d 

0 

-d 

, R = 

—r 

2 r 

—r 

0 

0 

0 

0 

0 

0 

—r 

r 

0 

0 

-d 

0 

d 


0 

0 

0 

0 


It turns out that for the pair ( D , R) condition (O fails to hold and hence by Theorem\l\ the oscillators 
do not synchronize. In particular, 


null 


R-rI q 

D 


1 ' 


' 1 ' 

0 

<£ range 

1 

-1 

1 

0 


1 


The solution corresponding to the eigenvector [10 — 1 0] T is shown in Fig. [3 where the first and third 
pendulums oscillate (with n radians of phase difference) at frequency ui = i/oJq +r while the second and 
fourth oscillators sit still. (We note that the two springs are identical.) Let us now figure out which 



Figure 3: Pendulums out of synchrony due to lack of observability. 


condition(s) of Assumption^ is violated for our example. First we consider the interconnection. The 
graphs Ts and Ta are given in Fig. [7} Since the graph Ts is connected, the third condition of Assumption]]) 
is satisfied. Note that Ta has two components: Ti and T 2 , the latter being a single vertex. Related to 
these graphs are the matrices R\ G R 3x3 and R 2 G R lxl . We have R 2 = 0 since T 2 has no edges. The 
matrix R± on the other hand has the following form 


Ri 


r —r 0 
—r 2 r —r 
0 — r r 
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Vi 




V3 



Ts 


ri 


r 2 






Figure 4: The graphs associated to the coupled harmonic oscillators in Example [T] 


For the graph T 2 the system da simply reads 1)2 + ^->oV 2 = 0 where r )2 = Z 4 £ R. Hence the as¬ 
sociated set of characteristic frequencies is singleton fl 2 = {wo}- The system associated to Fi reads 
rji + (ujqIs + Ri) rji = 0 wh ere r] 1 = [z\ Z 2 z^\ T £ R 3 . The set of eigenvalues of R\ being {0, r, 3r}, we 
have fli = {wo, v^o + r > \/ w o + 3r}. Now we can write fiinSl 2 = {wo}. Therefore the second condition 
of Assumption^ is also satisfied. Since the second and third conditions hold, the first condition must not 
(because the oscillators do not synchronize). The system ry 2 + WQ 772 = 0 is clearly observable from Z4. 
Therefore the other system iji + (wq I 3 + Ri)iji = 0 must be unobservable from at least one of its states 
Zk, k £ {1, 2, 3}. It can be shown that from Z 2 the system is indeed unobservable. This finding is not at 
all surprising when we look at the solution depicted in Fig. [3 


Remark 2 The situation shown in Fig. [3| not only renders the first condition of Assumption Q] more 
meaningful but also suggests a refinement on it. If the system in Fig. [1 were slightly modified by relocat¬ 
ing the damper between the third and fourth pendulums (as opposed to the original configuration where 
it connects the second and fourth pendulums) the observability condition of Assumption^ would still be 
violated yet the pendulums would this time synchronize. The reason is that even the component da 
described by the first three pendulums is unobservable from the second pendulum, it nevertheless is ob¬ 
servable from the third. And the significance of the third pendulum is that it is through it that the first 
component (in the modified system) is connected via damper to the second component (namely, to the 
fourth pendulum). One can carry this observation further so as to suggest the following relaxation of the 
first condition of Assumption |7J- “For all l = 1, 2 , ..., c the system (1131) is observable from each Zk for 
all k £ {i :vi £ Vi, d tj ^ 0 , j £ Ve}. ” 


Example 2 


Consider the following four coupled harmonic oscillators 
Z\ + ujqZi + r(zi - Z 2 ) = 0 

Z 2 + Wg^2 + d(z2 — £ 3 ) + r(z2 — Zi) = 0 

Z 3 + + d(z 3 - Z 2 ) + r(z 3 - Zi) = 0 

Z 4 + + r (^4 — 23 ) = 0 


where ojq, d, r > 0. The associated D and R matrices are provided below. 


' 0 

0 

0 

0 ' 


r 

—r 

0 

0 " 

0 

d 

-d 

0 

, R = 

—r 

r 

0 

0 

0 

—d 

d 

0 

0 

0 

r 

—r 

0 

0 

0 

0 


0 

0 

—r 

r 


It turns out that for the pair ( D , R) condition J3]) fails to hold and hence by Theorem [7] the oscillators 
do not synchronize. In particular, 


null 


R - 2 rl q 
D 


1 ■ 


' 1 ' 

-1 

(jL range 

1 

-1 

1 

1 


1 


The solution corresponding to the eigenvector [1 — 1 — 1 1] T is shown in Fig. [5| where the first and fourth 
pendulums make a synchronized pair and the second and third pendulums make another synchronized 
pair. These pairs oscillate (with ir radians of phase difference between pairs) at frequency ui = \/oJq +2r. 
(We note that the two springs are identical.) Let us now figure out which condition(s) of Assumption Q] 
is violated here. The graphs T s and T a are given in Fig. 0 Since the graph is connected, the third 
condition of Assumption^ is satisfied. Corresponding to the two components Ti andT 2 are the matrices 
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Figure 5: Pendulums out of synchrony due to extra common frequencies. 
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Figure 6: The graphs associated to the coupled harmonic oscillators in Example [2] 


i?i = i?2 



For the graph Fi the system (El reads iji + (wq/ 2 + Ri)r]i = 0 where rji = [zi Z 2 } T € K 2 . It can be 
shown that this system is observable from each Zk, k = 1, 2. Since i? 2 = R\, the same argument is valid 
also for the system El associated to r 2 . Therefore the first condition of Assumption]^ is also satisfied. 
This implies that the second condition cannot hold. Let us verify that the second condition does not hold. 
The set of eigenvalues of R\ being {0, 2 r}, we have fli = {wo, + 2r}. The equality i? 2 = R\ implies 
fl 2 = O-i. Hence n fl 2 = {wq, i/wq + 2r} / {w 0 } as expected. 


4 Electrical networks 


g 12 fe 



Figure 7: Three coupled LC oscillators. 


Hitherto the coupled pendulums were the workhorse in our study of the dynamics ©• There are 
however other important systems that share the same model; for instance, coupled electrical oscillators. 
Consider an array of identical LC oscillators where certain pairs (*, j) are coupled via an LTI resistor 
(with conductance gtj) or an LTI inductor (with inductance iij) or both; see Fig. [7] Writing Kirchhoff’s 
Current Law (KCL) at each node then yields the dynamics © in terms of electrical parameters: 

9 9 

cc 'zi + Zi T ^ ( 9ij(,Zi Zj) T ^ ^ hiji^Zi Zj ) — 0 , i 1) 2, ..., q (1*1) 

3 =1 3 =1 
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where Zj denote the node voltages, Cq > 0 and > 0 are respectively the capacitance and inductance of 
an individual oscillator, and hij = t~j when there is an inductor (i.e., tij 7 ^ 0 ) that connects the ith and 
jth nodes and hij = 0 otherwise. We work with passive components, i.e., gij = (jj t > 0 and hij = hji > 0. 
We take ga = 0 and ha = 0. Note that the interconnection of such an array can be represented by the 
admittance matrix 


Y(s) 


EjJ/i.i( s ) - 2 / 12 ( 5 ) 

-2/21 (s) Ej2/2j(s) 


-yiq(s) 
- 2 / 2 g (s) 


-y q i(s) - 2 / 92 ( 5 ) Ej2/9i( s ) . 


(15) 


where Uij(s) = gij + (£ijs) 1 = gij + hij/s is the admittance of the coupling between the ith and jth 
nodes. Define the symmetric positive semidefinite matrices 



E j 9ij 

—312 

-319 


E j hij 

—h \2 

— hlq 

G = 

-321 

E j 32 j • 

-329 

, H = 

—h 21 

Ej h 2i ■ 

— h - 2 q 



-392 

E j 9qj . 



~ hq 2 

E j hqj _ 


Note that Y(s) = G + s 1 H. Also note that G/cq and H/c 0 correspond to the matrices D and R of 
the array ©• In other words, G represents the dissipative coupling and H the restorative coupling. 
Given a matrix A £ C nxn let now Afc(A) denote the kth smallest eigenvalue of A with respect to the real 
part. That is, ReAi(A) < ReA 2 (A) < ••■ < ReA„(A). By Corollary |T| we can then state that the LC 
oscillators (fl4l) synchronize if A 2 (G) > 0, i.e., if the dissipative coupling graph is connected. Note that 
the condition A 2 (G) > 0 is only sufficient when H j 0. Now we point out an interesting extension of this 
inequality, which turns out to manifest itself in terms of the admittance matrix: 


Theorem 4 The LC oscillators (ED synchronize if and only if 


Re \ 2 (Y(Ju>)) > 0 for all to > 0 . 


(16) 


Proof. Without loss of generality take cq = 1. Then G and H correspond to the matrices D and R of 
the array ©■ Hence, by Theorem [1] the LC oscillators (TH1) synchronize if and only if 


null 


H-XI q 

G 


C range l q for all A € C . 


(17) 


Our task therefore reduces to establishing the equivalence of (fTCl) and (1171) . where Y(jco) = G+ (joe) -1 #. 
Let A £ C be an eigenvalue of Y(juj) and £ € C 9 be the corresponding unit eigenvector, i.e., Y(joj)f = A£ 
and ||£|| 2 = £*£ = 1, where £* is the conjugate transpose of £. We can write 

A = ?Y(ju>)£ 

= fiG-ju-'HK 
= CGf-jCHf/u. 

Since both G and H are symmetric positive semidefinite matrices we have Im A = < 0 and 

Re A = £*G£ > 0. Therefore no eigenvalue of Y(jui) can be on the open left half-plane. Also note that by 
construction G l q = 0 and Hl q = 0. Therefore Y(jco)l q = 0 and we can let Ai (Y(jw)) = 0 for all w > 0. 

Suppose now ED fails. This means that Re A 2 (Y(jui)) = 0 for some u > 0. There are two possibilities, 
one of which is: (i) Im A 2 (Y (jw)) = 0. In this case the eigenvalue at the origin is repeated and it must have 
at least two eigenvectors. (Otherwise l q would be the only eigenvector for the eigenvalue at the origin 
and there would exist a generalized eigenvector £ £ C q satisfying Y(joj)f = l q . But such £ could not 
exist because it would lead to the following contradiction: q = 1^ l q = l^V(jw)£ = (Y(joj) 1 9 ) T £ = 0.) 
Therefore we can find £2 ^ range l q satisfying T(jw )£2 = 0. This implies 0 = QY^w)^ = £^G £2 — 
j£|IT£ 2 /w. Since both G and H are symmetric positive semidefinite matrices we can deduce G £2 = 0 
and U£ 2 = 0. That is, null G fl null H D span {£ 2 }. Thus (1171) fails. Let us now consider the other 
possibility: (ii) ImA 2 (T(jw)) < 0. Then we can write A 2 (Y(juj)) = —j/3 2 for some #2 > 0. Let £2 £ C 9 
be the corresponding unit eigenvector, i.e., Y(ju)f 2 = and £^£2 = 1. Clearly, £ 2 ^ range l q . We 
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can write —j/3 2 = £ 2 ^C 7 W )£2 = ZvlG^ — jf 2 -^ 2 /^. This yields £ 2^2 = 0. Consequently, G£ 2 = 0 and 
—jj 3 2C2 = Y(ju)£ 2 = ~jHf2/0J. Therefore £2 has to be an eigenvector of H. In particular we can write 
null G l~l null (H — fawl q ) D span {£ 2 } and (1T71) once again fails. 

To show the other direction suppose this time that dm fails. Then we can find an eigenvector £2 ^ 
range l q that satisfies G£ 2 = 0 and H £2 = fa& for some fa £ K. Note that fa has to be real because it is 
an eigenvalue of H, a real symmetric matrix. Then we can write Y(ju>)£ 2 = (G— 2 = —jfato -1 ^- 
That is, A = —jfaui -1 is an eigenvalue of Y(joj). We also have Y(jcu)l q = 0. Therefore Re A 2 (Y(juj)) = 0 
and (flGl) fails. ■ 

Remark 3 Though a simple reexpression of Theorem [3 Theorem [^] is nevertheless significant (from the 
synchronization point of view) for it suggests a natural way of combining the two different interconnection 
graphs: the dissipative coupling graph and the restorative coupling graph. The result is a single graph with 
complex-weighted edges whose Laplacian is the admittance matrix Y(juj). 


2/31 (s) 



Figure 8: Three coupled electrical oscillators. 


We now briefly discuss how far the range of condition (fl6l) might extend. To this end we consider a 
generalization of the dynamics m employing the transfer function approach; see Fig. [8] In this general 
setting, each of q identical oscillator^ is assumed to consist entirely of resistors, capacitors, and inductors, 
all LTI and passive, i.e., with positive resistance, capacitance, and inductance values. We denote by yo(s) 
the admittance of an individual oscillator as seen from a given pair of terminals. One of those terminals 
is connected to the common ground and the other connects the oscillator to the rest of the network. In 
short, we represent each oscillator by an LTI passive one-port. As for coupling, the connection between 
a pair (i, j ) of oscillators is through also an LTI passive one-port (containing only resistors, capacitors, 
and inductors) with admittance yij{s). Note that yij(s) = Vji(s). As before, we take yu(s) = 0 and when 
there is no direct connection between the pair (*, j) we have yij(s) = 0. The overall interconnection gives 
us Y(s), the q-by-q admittance matrix (fl5l) . The array of coupled oscillators, considered as a whole, we 
denote by Af(yo{s), Y (s)). The network Af(y 0 (s),Y(s)) is said to synchronize if the node voltages z% (see 
Fig. @) synchronize, i.e., \zi(t) — Zj(t)\ —>• 0 for all i, j £ {1, 2, ..., q} and all initial conditions. (Note 
that here the initial condition, which determines the future evolution of the array, is the collection of all 
the initial capacitor voltages and initial inductor currents throughout the entire network.) In the sequel 
we will seek conditions guaranteeing the synchronization of A/"(yo(s), Y(s)). 

Thanks to passivity we will be able to proceed in our analysis fairly rapidly. First, note that the 
eigenvalues that are observable from the node voltages have to be the (finitely many) roots of the char¬ 
acteristic polynomial n(s) where n(s)/d(s) = det [yo(s)I q + Y(s)] and the polynomials n(s), d(s) are 
coprime. Since the network is passive those eigenvalues are confined to the closed left half-plane. Clearly, 
the ones with strictly negative real parts do not play any role in the steady state behavior of the network. 
This allows us to focus on the eigenvalues on the imaginary axis. Suppose now A = joj with uglis such 
an eigenvalue. Then (and only then) there exists a nonzero £ € C 9 (£ £ R 9 if lo = 0) and z{t) = Re(^e :,wt ) 
is a possible trajectory that can be traced by the node voltage vector z = [z± 22 • • • z q ] T . KCL imposes 

2 Although we stick to the term oscillator , the general setting we consider here allows also systems that do not display 
oscillatory behavior, e.g., a single capacitor. 
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on this trajectory the constraint Re ( [yo{jw)I q + Y= 0. In other words, 

£ e null [yo(ju)I q + Y (jw)] =: £(jui). 

Therefore any steady state solution can be written as a sum of finitely many terms z{t ) = Yh k R- e (£fc e:,li,,tt )> 
where oj k are distinct and £ £{juj k ) are nonzero. Evidently, this steady state solution corresponds to a 
synchronized collection of node voltages if and only if £ range l q for all k. Note also that £ range l q 
implies Y(ju k )f k = 0. Since £ k £ £(jw fc ), this means y 0 {j^k)kk = 0, i.e., y 0 {jio k ) = 0. Hence we 
obtained: 


Theorem 5 The network Af(yo(s), Y(s)) synchronizes if and only if 

null [yo(ju))I q + Y{j(jS)\ C range l q for all oj £ R. 

Also, for a synchronizing network, the steady state node voltages have the form Zi(t) = Y^, k ^ e ( a k^ Ukt ) 
with a k £ C and oj k £ R satisfying yo(jio k ) = 0. 


Purists may rightfully contend that the subspace £(jcv) is not always well-defined because for certain 
frequencies ui either yo(juj) or some entries yij(jai) of Y ( jui ) may attain infinite magnitude. This however 
is only a minor mathematical obstacle, easy to circumvent by thinking in terms of the physical system 
that the model stands for. Consider the case 12/o(J^) I = oo. This means that the impedance yff (jw) 
is zero, i.e., the oscillators behave as short circuit at that particular frequency ui. As a result, all nodes 
are grounded, i.e., all the node voltages Zi have to be zero, meaning £(Jlo) = {0}. Consider now the 
other potentially ambiguous case, \y%j{joj)\ = oo for certain pairs (*, j) while \yo{jui)\ < oo. Note that 
\yij(jaj)\ = oo implies that the nodes i and j are short-circuited, i.e., Zi = Zj. To get rid of the infinite 
terms yij(juj) = yjifjco) in the expression null [yo(ju))I q + Y(juj)\ we can remove ith and jth rows from 
the matrix [yo(jui)I q + Y(juj)\ and inject the following two new rows: (i) the sum of the removed pair of 
rows and (ii) a row that imposes the equality Zi = Zj. The cure can be repeated until all the infinite terms 
are gone. Let us demonstrate the procedure on an example network with q = 4 nodes whose admittance 
matrix reads 


Y{s) 


yi2{s) + 2/14 (s) 

- 2/12 (s) 

0 

-2/14 (s) 


— 2 /ia(s) 

2/12 (s) + 2/23 (s) 
“2/23 (s) 

0 


0 

-2/23 (s) 

2/23 (s) + 2/34 (s) 
“2/34 (s) 


—2/14(s) 

0 

2 / 34 ( 5 ) 

2/14 (s) +2/34 (s) 


Suppose that at some frequency w the admittances y 12 (ju>) and 2/23 (jw) are infinite. This gives us the 
equalities z\ = Z 2 and Z 2 = z 3 . Summing up the first three rows of \yo(joj)I q + Y{juj)\ lets us get rid of 
the terms 2/12 (jw) and 2 / 23 (jw). Then we inject the extra rows [1 — 1 0 0] and [0 1 — 1 0] to represent 
the relation z\ = Z 2 = 23 . Hence we can express £{jto) as 


£(ju>) = null 


2/0 (jw)+2/14 (jw) y 0 (jw) 
“2/14 (jv) 0 

1 -1 

0 1 


2/0(jw) + 2/34(jw) -2/14(jw) - 2/340‘w) 

-2/34 (jw) 2/0 (jw) + 2/14 (jw) + 2/34 (jw) 

0 0 

-1 0 


with righthand side cleansed of the infinite terms. Incidentally, another relevant point we want to make 
has to do with the eigenvalues of Y(joj). Since the roots of the polynomial p( A) = det [— XI q + Y(joj)\ 
are the eigenvalues of Y(ju>), we can use the procedure described above to define the finite eigenvalues of 
Y (jco ) when some of its entries are infinite. For instance, for the previous example, the finite eigenvalues 
of Y(ju> ), when 2/12 (jw) and 2/23 (jeo) are infinite, are defined as the roots of the polynomial 

-A + 2/340'w) - 2 / 14 (jw) - 2/34(jw) 

-2/34 (jw) -A + 2/14 (jw) + 2/34 (jw) 

0 0 

-1 0 


p(A) = det 


—A + 2/14 (jw) -A 

-2/14 (jw) 0 

1 -1 

0 1 


Being thus able to single out the finite eigenvalues allows us to continue to use the notation A k (Y(joj)), 
which will henceforth stand for the fcth smallest finite eigenvalue of Y(jco) with respect to the real part. 

Consider now a network Af(yo(s), T(s)) that does not synchronize. For this network Theorem [5] 
assures us that there exist a vector £ ^ range l q and a frequency ui £ K. satisfying Y{joj)f = —yo{j^)f, 
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and \yo{ju)\ < oo. Therefore A = —yo(i w ) is a (finite) eigenvalue of Y(ju>). Since the network is passive, 
all the eigenvalues of Y(juj) belong to the closed right half-plane, i.e., Re A k(Y(jco)) > 0 for all k, whence 
Re A > 0. Also, again due to passivity, Reyo(ju) > 0, whence Re A < 0. Consequently, Re A = 0. This 
implies, since Y(juj)l q = 0, the matrix Y(ju>) has at least two eigenvalues on the imaginary axis. This 
allows us to assert ReA 2 (Y(jw)) = 0. To summarize: 

Corollary 2 The network Af(yo(s), Y(s)) synchronizes if Re \ 2 (Y(jui)) > 0 for all oj £ R. 

5 Conclusion 

In this paper we studied the synchronization of identical (linear) pendulums coupled via dampers and 
springs. We first presented a necessary and sufficient condition for synchronization and then pointed 
out a sufficient set of conditions that may occasionally turn out to be easier to check than the former. 
Toward the end of the paper we applied the results obtained for pendulums to understanding better the 
collective behavior of coupled oscillators in LTI passive electrical networks. In particular, we established 
a relation between the second smallest eigenvalue of the node admittance matrix and the tendency of the 
individual systems to oscillate in unison. 
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